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Abstract 

Many simplified models, such as the celebrated Camassa-Holm equation, ad- 
mit peaked solitary waves. However, it is an open question whether or not such 
peaked solitary waves can be derived from the fully nonlinear wave equations. In 
this paper, a positive answer is given to this open question. A new approach is 
put forward to investigate the progressive waves with permanent form in finite 
water depth. Based on the symmetry, the new approach admits not only the 
traditional progressive waves with smooth crest, but also the solitary waves with 
peaked crest, and thus is more general. Using the so-called evanescent-mode 
base functions which decay exponentially in the horizontal direction, the peaked 
solitary surface waves are gained first by means of the linear wave equations, 
and then confirmed by the fully nonlinear wave equations. They have many un- 
usual, unique characteristics, which clearly indicate the novelty of them. Based 
on these peaked solitary surface waves, a few theoretical predictions and expla- 
nations to some natural phenomena are given. All of these might deepen our 
understandings and enrich our knowledge about solitary waves. 

Key words Solitary wave, peaked crest, progressive, fully nonlinear 

1 Introduction 



Since the solitary surface wave was discovered by John Scott Russell in 1834, various 
types of solitary waves have been found. The mainstream models of shallow water 
waves, such as the Boussinesq equation [1], the KdV equation [2], the BBM equation 
[3] and so on, admits dispersive smooth periodic and solitary waves of permanent 
form: the wave elevation is infinitely differentiable in the whole domain. Especially, 
the phase speed of these smooth water waves is closely related to the wave height: 
in general, a progressive wave with higher amplitude propagates faster than a lower 
one. Such kind of smooth periodic and solitary waves have been the mainstream of 
the teaching and investigating of water waves for quite a long time. 

However, in theory, the discontinuity of water wave elevation appears accidentally. 
It is well-known that the limiting gravity wave has a corner crest with 120 degree, as 
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pointed out by Stokes [1| in 1894. It is a pity that the importance of such kind of 
discontinuity seems to be neglected, mainly because Stokes limiting gravity wave |1] 
was regarded to hardly appear in practice. About one hundred years later, Camassa 
& Holm [5] proposed a model for the shallow water waves, called today the Camassa- 
Holm (CH) equation 

Ut + 2uJUx - Uxxt + 3mM^ = 2UxUxx + UUxxx, (1) 

where u{x,t) denotes the wave elevation, x and t are the spatial and temporal vari- 
ables, w is a constant related to the critical shallow water wave speed, respectively. 
As pointed out by Camassa & Holm [5], the CH equation ([T]) has the solitary wave 
when < < 1/2. Especially, when a; = 0, the CH equation ([1]) admits the peaked 
solitary wave in the closed-form 

u{x, t) = c exp(— |x — ct\), 

whose first derivative is discontinuous at the crest, where c denotes the phase speed. 
Unlike the KdV equation and Boussinesq equation, the CH equation ([1]) can model 
both phenomena of soliton interaction and wave breaking, as mentioned by Con- 
stantin [6]. Mathematically, the CH equation ([1]) is integrable and bi-Hamiltonian, 
thus possesses an infinite number of conservation laws in involution, as pointed out 
by Camassa & Holm [5]. In addition, it is associated with the geodesic flow on the 
infinite dimensional Hilbert manifold of diffeomorphisms of line, as mentioned by Con- 
stantin [6] . Thus, the CH equation ([1]) has many intriguing physical and mathematical 
properties. As pointed out by Fushsstainer [7], the CH equation ([1]) "has the potential 
to become the new master equation for shallow water wave theory" . 

In contrast to the above-mentioned theoretical results, the discontinuity widely 
appears in practical flows, such as dam break in hydrodynamics and shock wave in 
aerodynamics. In fact, such kind of discontinuous problems belong to the so-called 
Riemann problem [8l|9], a classic research field. Therefore, the discontinuity of wave 
elevation seems to be reasonable not only in mathematics but also in physics. 

Currently, the closed-form solutions of peaked solitary waves of the Boussinesq 
equation, the KdV equation, the BBM equation, and the modified KdV equation are 
found by Liao [TU]. Besides, it is also found by Liao [llj that the CH equation ([T]) 
admits peaked solitary waves even when u ^ 0. Besides, when u ^ 0, Kraenkel 
and Zenchuk [12] gave the explicit cusped solitary waves of the CH equation, called 
cuspon. The so-called cuspon is a kind of solitary wave with the 1st derivative going 
to infinity at crest. Note that the peakon has a finite 1st derivative, but the cuspon 
has an infinite 1st derivative. Thus, peakons and cuspons are completely different 
two kinds of discontinuous solutions of the CH equation ([1]). Therefore, nearly all 
mainstream models of the shallow water waves admit the peaked solitary waves. It 
indicates that the discontinuity and the peaked crest might be a common property of 
shallow water waves. 

Where does this kind of discontinuity come from? Are there any other peaked 
waves in finite water depth? It should be emphasized that all of the above-mentioned 
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mainstream models of shallow water waves are approximations of the fully nonlinear 
wave equations. Can we gain the peaked solitary waves from the exact nonlinear 
water wave equations? This is still an open question up to now. 

In this article, a positive answer to this open question is given. Using the fully 
nonlinear wave equations, we gain a new type of solitary surface waves in finite water 
depth, which have a peaked crest and possess many unusual characteristics quite 
different from the traditional ones. It is found that the fully nonlinear water wave 
equations admit two different kinds of waves: one is infinitely differentiable with phase 
speed closely related to the wave height, the other has a peaked crest. The former 
can be found in textbooks about water waves. However, the latter has never been 
reported, which suggests that the discontinuity might be a common property of water 
waves, and besides might enrich our knowledge and deepen our understandings about 
solitary waves. 

In §[21 the exact governing equation and boundary conditions for progressive waves 
with permanent form in finite water depth are described in a general frame, which 
admit not only the traditional periodic or solitary progressive waves with smooth crest 
but also the new type of peaked solitary waves. In § [3l a new type of peaked solitary 
surface waves are obtained by means of the linearized wave equations. In § |U the 
existence of such kind of new solitary surface waves with peaked crest is confirmed by 
the fully nonlinear wave equations. This kind of new peaked solitary surface waves 
have unusual characteristics different from traditional ones, as described in § 14.41 The 
concluding remarks, discussions and some theoretical predictions are given in § O 

2 Mathematical formulations 

First of all, we describe the mathematical formulations for progressive waves with 
permanent form in finite water depth in a general frame. Especially, we must be 
extremely careful so that the discontinuous solutions are not lost. 

Consider a progressive surface gravity wave propagating on a horizontal bottom 
with a constant phase speed c and a permanent form. For simplicity, let us consider 
the problem in the frame moving with the phase speed c. Let x, z denote the horizontal 
and vertical co-ordinates, with x = corresponding to the wave crest and z = —1 
corresponding to the bottom, i.e. the z axis is upward. Assume that the wave elevation 
has a symmetry about the crest at x = 0. Due to this symmetry, we need consider 
the domain x > only. Assume that the fluid is inviscid and incompressible, the 
flow is irrotational in the domain x > 0, and the surface tension is neglected. It 
should be emphasized that, due to the symmetry, the flow at x = is not necessarily 
irrotational. In other words, there might exist the vorticity at x = 0. In this way, a 
more generalized wave theory can be put forward, as shown below. 

Let D denote the water depth, g the acceleration due to gravity, the velocity 
potential (in the interval < x < -I-cxd), ( the free surface, respectively. All of these 
variables are dimensionless by means of D and y/gD as the characteristic scales of 
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length and velocity. Due to the symmetry C(~^) — C(^)) we only need consider the 
flow in the interval x e (0, +00), governed by 

VV(a;,^) = 0, ^ < C(a;),0 < x < +00, (2) 

subject to the boundary conditions on the unknown free surface z — ({x): 

«'^ + |^-«^ (V0-V0) + V0-vQv0-V0^ =0, 0<x<+oo, (3) 

C - + -V(/) • V0 = 0, 0<x< +00, (4) 
ox 2 



and the bottom condition 

d(f) 



= 0, z ^ -1, < X < +00, (5) 
oz 



where V is a Laplace operator and 



« = -F^ (6) 



is the dimensionless wave-speed, respectively Besides, we have either the periodic 
condition 

(f){x,z) = 4>{x + X,z), -l<z<C{x) (7) 

for the periodic progressive waves, where A = 27r//c is the wave length, or the decay 
condition 

(/)(±oo, ^) = 0, -l<z< C{x) (8) 

for the solitary waves. On the vertical boundary x = 0, we have the additional 
condition 

u{0,z) = hm^ = U{z), z<ax), (9) 

where U{z) is such an unknown horizontal velocity at a; = that the velocity po- 
tential 4>{x, z) and the corresponding progressive wave elevation C,{x) with permanent 
form exist. Furthermore, let Hy^ denote the dimensionless wave-elevation at x = 0, 
corresponding to the wave crest. For given i?^, one has an addition condition 

\imC{x) = H^. (10) 

In addition, the wave elevation must be bounded, i.e. 

|C(x)|<C, 0<x<+oo, (11) 

for a large enough constant C. Let u{x, y) and v{x, y) denote the velocity of the fluid. 
Due to the symmetry, in the domain x e (—00, -l-oo), we have the symmetry 

C{x) = C{-x), u{x, z) = u{-x, z), v{x, z) = -v{-x, z), (12) 
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which gives at x = that 

v{0,z) = ~v{0,z), i.e. v{0,z) = 0. (13) 

Since the flow is irrotational in the interval x G (0, +oo), the corresponding velocities 
u{x, z) and f (x, z) are given by 

. . d(f) . . d(f) 
u[x,z) = -—, v[x,z] = 7—, < X < +00. 
ox az 

Note that, at x = 0, we have f = due to the symmetry, and the boundary condition 
u = U{z) deflned by the limit ([9]) as x — )■ 0. It should be emphasized that, in order 
not to lose the solutions that are discontinuous at x = 0, the boundary condition IQ 
is deflned by a limit. In this way, the problem is deflned more generally than the 
traditional ones. 

Note that, according to the symmetry ( fT2l) about x = 0, the wave elevation ({x) 
and the horizontal velocity u are continuous at the vertical boundary x = 0. In fact, 
due to the symmetry (fT2|) . the boundary condition ([9]) is equivalent to the continuous 
condition of the horizontal velocity m at x = 0. It is a common knowledge that the 
Laplace equation ([2]) needs only one boundary condition at each boundary. Therefore, 
at X = 0, the boundary condition is sufficient for the Laplace equation ([2]) so that 
any other conditions for the smoothness of the horizontal velocity m as x — )■ are 
unnecessary: the higher-order derivatives of the horizontal velocity 

lim -— -, lim -— -, lim 



x^O dx"^ a;— >0 dx^ 2:— >-0 dx^ 

are unnecessary to be continuous, since the boundary condition ([9]) is enough for the 
Laplace equation ([2]). Thus, any other boundary conditions (such as that and C 
should be infinitely differentiable at x = 0) may lead to the loss of the discontinuous 
solutions and thus must be avoided. In other words, in the frame of the generalized 
water wave theory, both of (f) and C, are unnecessary to be infinitely differentiable at 
X = 0. 

The two nonlinear boundary conditions ([3]) and (jlj) must be satisfled on the 
unknown free surface z = C,{x). This leads to the mathematical difficulty to solve the 
nonlinear partial differential equations (PDEs). In case of small wave-amplitude, the 
linear boundary condition 

r.d'^S dd) , , 

a^T— ^ + -- = 0, onz = 0, 0<x< +cx), (14 
ox^ az 

is a good approximation of ([3]), and 



< X < +CX), (15) 



z=0 



is a good approximation of (jl]), respectively. The above two linearized free-surface 
boundary conditions, combined with the Laplace equation ([2]), the bottom condition 
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([5]) and either the periodic condition ([7]) or the decay condition ([8]), provide us the 
so-called linear wave equations. 

Based on the above linearized or fully nonlinear wave equations, hundreds of 
articles have been published for the periodic and solitary progressive waves. Nearly 
all of these traditional progressive waves are based on the base functions 

cosh[nA;(2; + 1)] sin(nA;x), n > 1, (16) 

for the velocity potential 0, which automatically satisfy the Laplace equation ([2]), the 
bottom condition ([5]) and the periodic condition ([7]), where k denotes the wave number 
and n > 1 is an integer. For periodic progressive waves with small wave-amplitude, 
substituting the velocity potential 

(j){x, z) = aAQCOsh[k{z + 1)] sin(A;x) (17) 

into the linear boundary condition f lT4|) . one has the dimensionless phase speed 



a=A/ ^<1' (18) 

say, the phase speed of a spatially periodic progressive wave in a finite water depth D 
is always less then \/ gD. Besides, substituting f|T7|) into f|T5|) gives the wave elevation 
with small amplitude 

C = ^ cos(fcx), (19) 

where = 2y4osinh(fc). The corresponding horizontal velocity reads 

aH^k cosh.\k(z + 1)] cos(kx) , , ,-, , , , 

u(x, z) = ^ \ ^ — = ^— - cosh A; z + 1 cos(kx), 20 

^ ' ^ 2sinh(A;) 2acosh(fc) ^ ^ ^ ^' ^ ^ 

which gives 

u cosh\k(z + 1)1 cos(A;a;) , , 

1T = JL ^ 0<x<+oo, 21 

Uo cosh(fc) 

where Uq = Hyj/{2a). At x = 0, we have the corresponding horizontal velocity 

H U 

Uiz) = lim u(x, z) = ^ cosh[k(z + 1)1 = — tttt cosh[k(z + 1)1 . 

x^o 2acosh(/c) cosh(/c) 

Note that the velocity potential and the wave elevation given by the above traditional 
linear wave theory automatically satisfies the symmetry f[T^ . In other words, we 
gain exactly the same results by first solving the PDEs to flTT]) in the interval 
< X < +00 and then expanding the result to the interval — oo < x < by means of 
the symmetry (fT2l) . At x = 0, the horizontal velocity u = U{z) is gained by the limit 
([9]), and the vertical velocity is continuous, since 



(22) 
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which is equal to the boundary condition v = that is based on the symmetry (fT2|) . 
Therefore, the mathematical formulas ([2]) to (fTTl) are consistent with the traditional 
wave theory, although it is more general. 

Note that, for given x, the horizontal velocity u of periodic progressive Airy's 
waves decreases exponentially as z varies from the surface {z = 0) to the bottom 
{z = —1). Especially, based on the base functions (1T6|) . the elevation of the Airy's 
wave and the corresponding velocities are infinitely different iable, although such kind 
of smoothness is unnecessary at all for the generalized wave theory described by ([2]) 
to (HI]). 

For the periodic progressive surface waves with large amplitude, the fully nonlin- 
ear wave equations must be considered. As pointed out by Cokelet [13], the phase 
speed c of the progressive periodic waves depends not only on the water depth D and 
the wave number k but also on the wave height H^'- in most cases, the larger the wave 
amplitude, the faster the periodic wave propagates. In other words, the traditional 
progressive periodic waves are dispersive. Besides, the periodic progressive surface 
waves have a smooth crest with the exponentially decaying velocity u{x, z) from the 
surface to the bottom. Like the Airy's linear waves, the traditional nonlinear periodic 
progressive waves are also infinitely different iable, although boundary conditions for 
such kind of smoothness does not exist at all. Here, it should be emphasized that, the 
same results for the traditional progressive periodic waves can be obtained by first 
solving the PDEs (|5]) to f lTT]) in the interval < a; < +oo and then expanding the 
results to the interval —oo < a; < by means of the symmetry f[T^ . Besides, at x = 0, 
the horizontal velocity u given by (Q and the vertical velocity f = agree completely 
with the traditional ones. Thus, Eqs. ([2]) to ( ITT]) indeed admit the traditional waves 
with smooth crest. 

It should be emphasized that, in the frame of the linear wave theory, solitary 
waves have never he reported, to the best knowledge of the author. For details, please 
refer to Mei et al . Solitary wave solutions for nonlinear and dispersive long waves 
had been found by Boussinesq [1] and Rayleigh |15]. For dispersive long waves of 
permanent form, the so-called KdV equation [2] gives the periodic cnoidal wave for a 
finite wavelength A, which tends to the solitary wave 



((x) = Hu, sech^ 

with the phase speed 

c = + H^, (24) 

as A — i- +00. Note that these cnoidal and solitary waves have a smooth elevation, 
say, C,{x) is infinitely differentiable for all x G (— oo,+oo). Besides, its phase speed c 
depends upon the wave height Hyj-. the larger the wave height, the faster the solitary 
wave propagates, as shown by fl24]) . All of these results can be gained by first solving 
the KdV equation in the domain < x < +oo and then expanding the results to the 
domain — oo < x < by means of the symmetry ( |T2i) . 

Both of the above-mentioned solitary wave and the Airy linear waves are special 



[X 



ct) 



(23) 
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cases of the so-called cnoidal waves. By means of perturbation methods and using the 
fully nonlinear wave equations, Fenton [T6l[T7] gave respectively a high-order cnoidal 
wave theory and a ninth-order solution for the solitary wave in the form 



where ajj,e,/3 are constants determined by the physical parameters. It should be 
emphasized that all of these traditional cnoidal and solitary waves have a smooth 
crest: ({x) is infinitely differentiable for all x G (—00, +00). Besides, the velocity 
u{x, z) at bottom is always larger than that at crest. Furthermore, the phase speed is 
dependent upon wave height. Finally, to the best of author's knowledge, all traditional 
solitary surface waves have a crest higher than the still water: the solitary waves in 
the form of depression have been reported for interfacial waves, but never for the 
surface waves. It should be emphasized that all of these traditional results can be 
gained by first similarly solving the PDEs ([2]) to (|TT1) in the domain < a; < +00 and 
then expanding the results to the domain —00 < x < by means of the symmetry 
f[T^ . All of these indicate that the PDEs ([2]) to flTT]) with the symmetry condition 
(fT2|) are indeed consistent with the traditional linear and nonlinear progressive waves 
with smooth wave crest. 

Indeed, the above-mentioned traditional periodic and solitary progressive waves 
are infinitely differentiable. However, this kind of smoothness is unnecessary, since no 
such kind of smoothness conditions are enforced to the PDEs ([2]) to flTT]) . In essence, 
such kind of perfect smoothness of the wave elevation and velocities automatically 
come from the base functions (fT6|) . which are infinitely differentiable at x = 0. 

There exist a little thing in the traditional theories of solitary wave that should 
be reconsidered carefully. Note that the traditional cnoidal waves are periodic and 
thus have an infinite number of wave crests. As a special case of the cnoidal waves 
as the wave-length A — )■ +00, the solitary waves f l2^ of the KdV equation should 
have an infinite number of wave crest, although the distance between the two crests 
is infinite. So, seriously speaking, the solitary waves ( 123|) given by the KdV equation 
is not truly "solitary" , since it might have an infinite number of crests. 

Note that, like the base functions f fT6|) that are widely used for the smooth periodic 
and solitary progressive waves, the following base functions 



automatically satisfy the Laplace equation ([2]), the bottom condition ^ and the 
bounded condition f llip . However, different from the smooth base functions f ll6p . the 
above base function decays exponentially in the x direction and satisfies the boundary 
condition (|8]). Thus, it is more convenient to strictly express solitary waves that have 
truly only one crest. Especially, unlike the base function f|T6|) . its derivatives with 
respect to x are not differentiable at a; = 0. 

Even so, the base function with discontinuity of the 1st derivative at crest was 
widely used as the so-called evanescent (or non-propagating) mode [TS] in the problems 




(25) 



i=i j=i 



cos[nA;(z + 1)] exp(— nA;x), n>l, A;>0, 0<x< +00 



(26) 
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of linear waXer wave diffraction-refraction by discontinuous bed undulations [T91422] . 
or linear waves propagating over a bed consisting of substantial variations in wa- 
ter depth p3]425] . and so on. The solutions of these problems contain not only the 
propagating waves expressed by the base function ( IT6l) . but also the non-propagating 
(or evanescent) waves expressed by the base function (l26l) . which represent localized 
effects and depend on the local bottom geometry [25]. However, to the best of the au- 
thor's knowledge, the base function f l26|) has never been applied to express progressive 
solitary waves with permanent form propagating in a constant water depth. 

As mentioned above, mathematically speaking, no smoothness conditions are en- 
forced to the PDEs ([2]) to (fTTl) . Physically, such kind of discontinuity widely appears 
in practice, such as dam break and shock waves, which have clear physical meanings. 
Thus, like the peaked solitary waves of the CH equation ([1]), the peaked solitary waves 
of the PDEs ([2]) to f lTT]) should be reasonable and acceptable not only in mathematics 
but also in physics. 

Can we find any kinds of solutions of peaked solitary surface waves of the fully 
nonlinear wave equations ([2]) to (ITT]) by means of the evanescent base function (12^ ? 
The answer is positive: we demonstrate in the following part of this article that the 
fully nonlinear wave equations 02]) to f lTT]) indeed admit such a new type of solitary 
surface waves with peaked crest and some unusual characteristics that are completely 
different from the traditional ones. 



3 Solitary waves by linear equations 

As mentioned in § 2, both of the base functions f|T6]) and fl26]) automatically satisfy 
the Laplace equation ([2]), the bottom condition ([5]) and the bounded condition f lTT]) . 
In the frame of the linear wave theory, the former satisfies the periodic boundary 
condition (JT]) and gives the well-known Airy wave, which is infinitely differentiable and 
horizontally periodic, as mentioned above. The latter satisfies the decay condition (JH]) 
and gives peaked solitary waves, as shown below. 

In the interval < a; < +oo, we have the velocity potential in the form 

(f)+{x,z) = aAcos[k{z + 1)] e''''', < x < +oo, (27) 

where the superscript + denotes the interval x G (0, +oo), a denotes the phase speed, 
A is a constant related to the wave height, and > is a parameter determined by the 
phase speed a, respectively. Note that the above expression automatically satisfies 
the Laplace equation the bottom condition the decay condition and the 
bounded condition (fTT]) . Substituting ( 127]) into the linear boundary condition (IT^ 
gives 

akA (a^/c cos k — sin A;) exp(—kx) = 0, < x < +oo, 

which leads to 

9 tank , TT , , 

a = , niT < k <n7f + -, (28) 
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where n > is an integer. Similarly, defining 0~ in the interval x G (—00, 0), we gain 
exactly the same result. Obviously, a > 1, i.e. c > y/gD. This is different from the 
periodic Airy wave whose phase speed has the property c < y/gD. 

Given the dimensionless phase velocity a, the above equation has an infinite 
number of solutions /c„, where 

nvr < A;„ < nvr + ^, n > 0. (29) 

For the sake of simplicity, define the set 

Kq, = <kn ■ = — -J — nvr < /c < nTT + — , n = 0, 1, 2, 3, ■ ■ • > . (30) 

Write kiy G Kq, where i/ > is an integer. Then, we have the solution 

4>^{x, z) = aAy cos [ky{z + 1)] exp(— A;,^x), < x < +00. (31) 



According to the linearized boundary condition f[T5|) . the corresponding elevation of 
the solitary wave reads 



(+(x) = a 



= —a A^k„ cos k,y exp{—k„x) 

z=0 



dx 

= —A^smk^exp{—k,^x), < x < +cxd, (32) 

where a'^k^, cos k^, = sin k^, due to fl28l) is applied. Due to the continuity of wave 
elevation, we have at x = the wave elevation 

C(0) = limC+(x). (33) 

Then, according to (fTOl) . we have the relation 

H^ = —A^sink^. (34) 

Thus, using the symmetry f|T2|) . we have the peaked solitary wave 

C(x) = i^^ e"^'"l''l, -cx) < X < +00. (35) 

This is a solitary wave that seriously has only one crest, with a discontinuous 1st- 
derivative C'i^) crest! For example, in the case of a = 3^/^/i/7r, we have ko = 
7r/3, ki = 4.58117, and the corresponding peaked wave elevations are shown in Fig. [H 
respectively. All of them have a peaked crest. Besides, the larger the value of k,y, the 
sharper the peaked solitary wave. These are essentially different from the traditional 
periodic and solitary progressive waves, which are infinitely differentiable. This clearly 
indicates the novelty of the new type of peaked solitary waves. 

In the interval < x < +00, the corresponding horizontal velocity reads 

90+ ak^H^cos[k^{z + l)]e-''-''' 

u+{x,z) = -7^ = r^T^ . 36 

ox sm{ki,) 
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Figure 1: Q{x)/H^ in the case of /cq = 7r/3,/ci = 4.58117 with a = S^/'^/a/tt- 
Dashed hne: exp{—kQ\x\); Sohd hne: exp{—ki\x\); Dash-dotted hne: [exp(— A;o|a;|) + 
exp{—ki\x\)]/2. 




Figure 2: Velocity profile u/Uo at x = in the case of ko = n/S with Uq = H^/a. 
Solid line: periodic Airy wave; Dashed line: linear peaked solitary wave. 
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Using the symmetry (fT2|) . we have 

u [x,z) = u^{—x,z) = — , — oo < a; < 0. (37) 

sin(fc,,) 

At X = 0, M is continuous and we gain the corresponding horizontal velocity 

U[z) = hm — — = . . . (38) 

Thus, in the whole domain — oo < x < +oo, we have a uniform expression 

u cos[fc,,(z + l)]e"''"l^l 



Uq cos{ki, 



X G (— cxD, +oo), (39) 



where Uq = H^/a. So, for given a;, the horizontal velocity u of the peaked solitary 
wave increases as z varies from the surface {z = 0) to the bottom [z = —1): in 
other words, u at bottom is always greater than that on surface. For example, when 
ko = 7r/3 (corresponding to the phase speed c = 3'^^'^ / y/ii), the horizontal velocity at 
bottom beneath crest of the peaked solitary wave is twice of that on surface, as shown 
in Fig. [21 This is quite different from the traditional ones whose horizontal velocity 
u at bottom is always less than that on surface. This also indicates the novelty of the 
new solitary waves. 

Let v~^{x,z) and v~{x,z) denote the vertical velocity in the interval x > and 
X < 0, respectively. In the domain < x < +oo, the corresponding vertical velocity 
reads 

A - ^ - «fe^-g^sin[fc,(z + l)]e-^-" 
^ ~ dz ~ sinfc, ■ ^ 

Using the symmetry (fT2|) . we have 

_ +1 \ akyH^sin[ku{z + l)]e^'''' 

V (x, z) = —v^[—x,z) = — ; , — oo < X < 0. (41) 

sin ku 

As X — )■ 0, we have the limits 

lim = ^ sm.[ky{z + 1)1, hm v~ = — ^^^ ^ sm.[ky{z + 1)1 (42) 

x-j-o sin ky sin k^ 

so that 

lim = — lim v~ ^ 0. (43) 

However, according to f|T3l) . the vertical velocity v exactly equals to zero at x = 0, 
i.e. f (0, z) = 0. Thus, at the interface x = 0, the vertical velocity v is discontinuous: 
V changes sign as we cross the interface x = 0. Besides, for given z, the jump of v at 
X = 0, i.e. 

lim(f ^ — v^) = 2 lim , 
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is directly proportional to kiy. It should be emphasized that, in the frame of the 
traditional wave theory, the discontinuity and jump of velocity on an interface is 
familiar and acceptable. For instance, on the free surface of interfacial waves [26H29], 
although the velocity normal to the interface is continuous, the tangential velocity is 
discontinuous, as pointed out by Lamb [30] (§231, page 371): "the tangential velocity 
changes sign as we cross the surface" , and "in reality the discontinuity, if it could ever 
be originated, would be immediately abolished by viscosity" . So, the discontinuity of 
f at X = in the frame of the generalized wave theory is consistent with the traditional 
wave theory. 

It should be emphasized that, given a dimensionless phase velocity a, there exist 
an infinite number of ki, satisfying (1281) , and each of them corresponds to a solution of 
the linear wave equations. This is quite different from the smooth periodic progressive 
Airy wave ({x) = Acos{kx) that is unique for a given a. This is because, given a 
dimensionless phase speed a, the transcendental equation f|T8l) has an unique solution, 
but the transcendental equation (1251) has an infinite number of solutions! Note that, 
unlike the smooth Airy wave theory, the parameter k^, of the peaked solitary wave 
does not denote the wave number, but the decaying rate of the wave elevation as 
X —7- +oo: the larger the value of k^, more quickly the wave elevation decays to zero. 
Besides, according to (128|) . the peaked solitary wave can propagate very quickly even 
if the wave height is small, since tan. k^y/k^, — > +oo as k^y — )• z/vr + 7r/2, where 
ki, G Kq and z/ > is an integer. 

Since the governing equation ([2]), the linearized free surface conditions f[T^ and 
(fT5|) . and all other conditions are linear, the solution of the peaked solitary wave can 
be expressed in a general form 



Using the linear boundary condition (ITBI) . we have the corresponding elevation of the 
solitary wave 



+ 00 

(p~^{x, z) = a An COS [kn{z + 1)] exp(— /c^a;), Q < x < +oo. 

n=0 



(44) 



a 



(90 



dx 



,+ 



= —a 



+ 00 

^ ^ Ankn COS kn exp( knX) 
n=0 




n=0 




(45) 



n=0 
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where the relationship (!28|) is apphed. Using the symmetry (fT2!) . we have 

+00 

Ci.^) = ^bnexp{-kn\x\), (46) 

n=0 

+00 

u{x,z) = -a'^AnknCos[kniz + l)]exp{-kn\x\), (47) 

n=0 

+00 

-a ^ A„ fc„ sm[kn{z + 1)] exp(-A;„a;), (48) 

n=0 

+00 

: ^ A„ sin[A;„(z + 1)] exp(A;„x), (49) 



V [X,Z) 



V [x,z) = a 



n=0 



where the superscripts + and — denote the interval x G (0, +00) and x G (— oo,0), 
respectively. Note that, using the restriction ffTOj) . we have a linear algebraic equation 



+00 



^«, = 2^fen, (50) 
n=0 

which has an infinite number of solutions of bn for a given value of Hy^, where n > 
is an integer. In other words, given a wave height and a phase speed c, there 
exist an infinite number of peaked solitary waves. This is quite different from the 
smooth periodic progressive Airy wave that is unique for given phase speed and wave 
height. This is because, given a phase speed a, the transcendental equation ( IT8l) has 
an unique solution but ( 128|1 has an infinite number of ones! Thus, different peaked 
solitary waves may have exactly the same phase speed. For example, in the case of 
ko = 7r/3, ki = 4.58117 when a = S^/^/yvr, the peaked solitary wave 

is different from exp(— fcol^^l) and exp(— A;i|x|) but has the same phase speed with 
them, as shown in Fig [1] Note that, in the frame of the generalized linear wave 
theory, the crest of each peaked solitary wave may be at different position. So, the 
peaked wave elevation can be generally expressed by 

+00 

Cix) = ^bnexp{-kn\x - ^o\), (51) 

n=0 

where 6„ and ^„ are real numbers. 

Like Airy's wave, since the elevation fl5S]) of the peaked solitary wave is gained 
by the linear wave equations, the value of can be negative, corresponding to a 
peaked solitary wave in the form of depression. For example, ({x) = — exp(— |x|)/10 
is a peaked solitary wave of depression. Such kind of peaked solitary waves have been 
never reported for surface waves. This once again indicates the novelty of the new 
solitary waves. 



15 



It is traditionally believed that solitary waves are always governed by nonlinear 
differential equations. However, we illustrate here, for the first time, that the solitary 
waves exist even in the frame of the linear wave equations in finite water depth! Note 
also that the peaked solitary wave (l35ll is the same as the peaked solitary wave found 
by Casamma & Holm [5]. This might reveal the origin of the peaked solitary waves of 
the CH equation ([1]), since it is an approximation of the fully nonlinear wave equations 
in shallow water. However, the peaked solitary wave fl35|) is valid not only in shallow 
water but also in finite water depth, and thus is more general. 

Finally, it should be emphasized once again that these exists the discontinuity of 
the vertical velocity at x = 0, say, 

f (0, z) = but lim = — lim v~ ^ 0. 

Thus, V changes sign as we cross the interface a; = 0. However, such kind of dis- 
continuity is familiar and consistent with traditional wave theories. For example, as 
mentioned by Lamb [SU] (§231, page 371), "the tangential velocity changes sign as 
we cross the surface" of interfacial waves, and "in reality the discontinuity, if it could 
ever be originated, would be immediately abolished by viscosity" . It seems that such 
kind of velocity discontinuity is inherent in the frame of the new wave theory based 
on the evanescent-mode base-function ( I26p . which is however completely consistent 
with the traditional wave theory based on the smooth base-functions (fTB]l . 

4 Solitary waves by nonlinear equations 

As shown in § [31 the peaked solitary surface waves given by the linearized free surface 
conditions f[T^ and f lT^ has some unusual characteristics quite different from the 
traditional periodic and solitary ones. Does the fully nonlinear wave equations (|2]) 
to ( ITTj) indeed admit such kind of peaked solitary waves? Does this kind of peaked 
solitary waves have the same unusual characteristics as those given by the linear wave 
equations, if the answer of the above question is positive? 

To answer these questions, we consider here the solitary surface waves with finite 
wave-amplitude so that the nonlinear terms of the free surface conditions and 
are not negligible, and besides 2; = is not a good approximation of the free surface 
({x). In other words, we had to solve the fully nonlinear wave equations (j2]) to (fTTj) 
accurately. 

In this paper, an analytic technique for highly nonlinear differential equations, 
namely the homotopy analysis method (HAM) proposed by Liao [^TH5B] . is applied 
to solve the fully nonlinear wave equations to flTT]) . Unlike perturbation techniques, 
the HAM does not need any assumptions of small physical parameters, since it is based 
on the homotopy, a basic concept in topology. Besides, the HAM provides us great 
freedom to choose base functions for solutions of considered nonlinear equations. Espe- 
cially, by means of the so-called "convergence-control parameter" that has no physical 
meanings, the HAM provides us a convenient way to guarantee the convergence of 



16 



approximation series: in essence, it is the so-called "convergence-control parameter" 
that differs the HAM from all other analytic approximation techniques, as pointed 
out currently by Liao [36]. Therefore, the HAM is valid for highly nonlinear problems, 
as shown by many successfully applications in fluid mechanics, applied mathematics, 
physics and finance. For example, by means of the HAM, Liao [32] gained, for the 
first time, convergent series solution for Blasius and Falker-Skan boundary-layer flows, 
which are uniformly valid in the whole field of flow (Note that the traditional power 
series given by Blasius [37J is valid only in the near field, and thus had to be matched 
with another asymptotic approximation of flow in far field). Besides, using the HAM 
as a tool, the exact Navier-Stokes equations were solved by Turkyilmazoglu |3B] for a 
compressible boundary layer flow due to a porous rotating disk, and by Xu et al [39] 
for viscous flows in a porous channel with orthogonally moving walls. Furthermore, 
the limit cycle of Duffing - van der Pol equation was solved by Turkyilmazoglu [lO] , 
and the two coupled Van der Pol equations were solved by Li et al [4]J . Especially, by 
means of the HAM, some new boundary layer flows have been found by Liao [12] and 
by Liao & Magyari [13] , which have been neglected by other analytic and even numer- 
ical techniques. In addition, the HAM has been also successfully applied to solve some 
nonlinear PDEs with moving boundary conditions, such as those about American put 
option. For example, Zhu [44j successfully applied the HAM to give a series approxi- 
mation of the American put option, which gives optimal exercise boundary valid for 
a couple of years, while perturbative and/or asymptotic formulas are accurate only in 
a few days or weeks. All of these illustrate the potential and validity of the HAM for 
highly nonlinear problems. 

It should be emphasized that the HAM has been successfully applied to solve 
the fully nonlinear wave equations. Using the traditional base functions flTBl) . Liao & 
Cheung ^45J applied the HAM to solve the periodic progressive surface waves in deep 
water and obtained convergent solutions for waves with high amplitude even close to 
the limiting case. Their analytic results agree quite well with those given by Schwartz 
[i6] and Longuet-Higgins [U]. Besides, using the same traditional base functions (fT6|) . 
Tao et al [18] successfully applied the HAM to solve the fully nonlinear wave equations 
for periodic progressive waves in finite water depth, and their analytic results agree 
well not only with the analytic ones given by Cokelet [13] and Fenton ^S] but also 
with the experimental ones reported by Mehaute et al [50]. Currently, Xu et al [51] 
successfully applied the HAM to investigate the steady-state fully-resonant progressive 
waves in finite water depth, and found, for the first time, the multiple steady-state 
fully resonant waves without exchange of wave energy. All of these demonstrate the 
validity of the HAM for the fully nonlinear wave equations ([2]) - ([TTj) . 

4.1 Analytic approach based on the homotopy analysis 

As shown below, the generalized, fully nonlinear wave equations ([2]) to (ITT!) can be 
solved by means of the HAM and the evanescent-mode base functions (1261) in a similar 
way as those by Liao & Cheung [IS] and Tao et al [IS], although they used the 
traditional base functions f|T6l) . 
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Due to the symmetry (fT2|) . we need consider the case x > only. As mentioned 
in §[3] for the hnearized wave theory, given a phase speed a, there exists an infinite 
number of peaked sohtary waves ({x) = Hyjexp^—k^x) with different decay-rate ki, 
satisfying the transcendental equation (l28l) . where k^ > is an integer. So, the 
general expression of the peaked solitary waves becomes rather complicated when 
the nonlinear boundary conditions are considered. Thus, for the sake of simplicity, 
we consider in this section the peaked solitary waves only with the primary decay- 
rate < fco < vr/2. Here, let k denote the primary decay-rate fco, if not mentioned. 
Obviously, due to the nonlinearity of the two free surface conditions, the peaked 
solitary wave elevation should contain the terms exp(— n/cx), and correspondingly the 
velocity potential function (j){x, z) should be expressed in the form 

+00 

4){x,z) = ctn cos[nk{z + 1)] exp{—nkx), x > 0, k > 0, (52) 

n=l 

which automatically satisfies the governing equation ([2]), the bottom boundary con- 
dition ([5]) and the bounded condition ffTTj) . where k = ko E is defined by fl30|) . and 
(Xf^ IS a constant to be determined. We search for the solitary surface waves in the 
form 

+ CX1 

C(x) = bn exp {—nkx) , a; > 0, /c = /cq ^ K^,, (53) 

n=l 

where hn is a constant coefficient to be determined. The above expressions (!52l) and 
(153|) provide us the so-called solution-expression of 0(x, z) and respectively, 
which play important role in the frame of the HAM, as shown below. 

Let (poiyX , z) denote the initial guess of the velocity potential (/>(x, z) and 
the wave elevation C{x) in the interval x G (0, -|-oo), respectively. To apply the 
HAM, we should first of all construct two continuous variations from the initial guess 
0o(a^, z), Co{x) to the exact solution 0(x, z), C(a;), respectively. This can be easily done 
by means of the homotopy, a basic concept in topology, as shown below. 

First, according to the solution expression fl52|) . we choose 

0o(x, z) = Ao cos[k{z + 1)] e"^"", x > 0, A; > 0, (54) 

as the initial guess of the velocity potential 0(x,^), where Aq is a constant to be 
determined later. Note that, different from Liao & Cheung [35] and Tao et al [48], the 
evanescent-mode base function fl26|) is used here. Note also that (j)o{x, z) automatically 
satisfies the Laplace equation ([2]), the bottom condition ^ and the bounded condition 
ffTT]) . Besides, following Liao & Cheung [IS] and Tao et al [IB], we choose 

Co(x) = (55) 
as the initial guess of wave elevation C,{^)- 

Secondly, according to ([3]), we define a nonlinear operator 

86 d /I \ 

•^'^ = + ^ - «^ (V0 ■ V0) + V0 ■ V f -V0 ■ V0 j . (56) 
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Let q G [0, 1] denote an embedding parameter, c,^ and c,, be two non-zero auxiliary 
parameters without physical meanings, called the convergence-control parameters, 
and C denote an auxiliary linear operator, respectively. Following Liao & Cheung [45] 
and Tao et al [IE], we construct the so-called zeroth-order deformation equation 



V^$(a:, z; g) = 0, < s < +oo, z<r]{x;q), (57) 
subject to the boundary conditions on the unknown free surface z = r]{x; q): 

{1 - q)C[^{x,z;q) - (f)o{x,z)] = q M[^{x,z;q)] , < x < +oo, (58) 



(1 - q)7]ix;q) =c^q 
and the boundary condition at the bottom 



■ , 9$ 1 
r]{x] q) - a— + -V$ ■ V$ 
ox 2 



<x < +00, (59) 



9$ 

— = 0, z = -1, < X < +00. (60) 

OZ 

If wave height is given, there exists the additional condition: 

\imr]{x;q) = H^. (61) 

Note that $(x, z;g) and ri{x;q) depend not only on the original physical variables 
X, z but also on the embedding parameter g G [0,1] and the two convergence-control 
parameter c^, that have no physical meanings at all. It should be emphasized that 
we have great freedom to choose the values of the convergence-control parameters 
and Cj^. Following Liao & Cheung |15] and Tao et al [IE], we choose the auxiliary 
linear operator 

which has the property C[0] = 0. Note that C is exactly the linear part of the nonlinear 
operator defined by (156|) . In this way, the zeroth-order deformation equations ([5j 
to fl6T]) are well defined. 



When g = 0, we have from ( l59l) that 

r]{x;0)=0 = Coix), (63) 
and then the corresponding zeroth-order deformation equations become 

V^<^{x,z;0) = 0, z<0, 0<x<+oo, (64) 
subject to the boundary conditions on the known free surface 

C[^{x, z;0) - (f)o{x, z)] = 0, when 2; = 0, < x < +00, (65) 
and the boundary condition at the bottom 

^%^ = 0, z = -l, 0<x<+oo. (66) 

OZ 
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Since the auxiliary linear operator L has the property £[0] = and besides the 
initial guess 0o(a^, A defined by (15^ satisfies the Laplace equation ([2]) and the bottom 
condition ([5]), it is straightforward that 

<l>(x,2;;0) = </)o(a;,z). (67) 

When g = 1, since ^ ^ and Cr^ ^ 0, the zeroth-order deformation equations 
( l57|) to ( !6T]) are equivalent to the original fully nonlinear wave equations ([2]) to (fTOl) . 
respectively, so that we have the relationship 

^{x,z-l) = (t){x,z), r/(a;;l) = C(a;). (68) 

Thus, as the embedding parameter q increases from to 1, $(a;, z;g) and ri{x;q) 
indeed vary continuously from the initial guess (f)o{x , z) , (oi^) to the exact solution 
0(x, z), ({x) of the fully nonlinear wave equations ([2]) to (ITU]) , respectively. Therefore, 
the zeroth-order deformation equations fl57|l to fl6T|) truly construct such a kind of 
continuous variation that provides a base of our analytic approach, as shown below. 

Since both of $(x, z; q) and ri{x; q) are dependent upon the embedding parameter 
q G [0,1], we can expand them in Maclaurin series with respect to q to gain the 
so-called homotopy-Maclaurin series 

+ 00 

^{x,Z]q) = Mx,z) + ^(l)m{x,z) q"", (69) 

m=l 

+ 00 

vix;q) = $^C™(x)g^ (70) 

1 d'^r](x; q) 



m=l 



where 



1 9-$(x,2:;g) 

(pm{x,Z) - 



ml dq^' 



Cm('^) 



ml dq^' 



g=0 



g=0 

and the relationship ( l63ll and (1671) are used. However, it is well known that a Maclaurin 
series often has a finite radius of convergence. Fortunately, both of $(0;, 2;; q) and 
ri{x\ q) contain the two convergence-control parameters and c,,, which have great 
infiuence on the convergence of the Maclaurin series of $(a;, z; q) and 77(0;; g), as shown 
by Liao & Cheung and Tao et al [48j. Here, it should be emphasized once again 
that we have great freedom to choose the values of and c^. Thus, if the convergence- 
control parameters c^, are properly chosen so that the above homotopy-Maclaurin 
series are convergent at g = 1, we have the homotopy-series solution 

+00 

0(X,Z) = (t)o{x,z) + ^(t)m{x,z), (71) 

m=l 

+00 

C(x) = ^C™(x). (72) 



m=l 



The equations for the unknown (j)m{x,z) and Cm(a^) can be derived directly from 
the zeroth-order deformation equations. Like Liao & Cheung flS] and Tao et al [48], 



20 



substituting the series (169|1 and (170!) into the zeroth-order deformation equations (!57|) 
to (16T|) . then equating the hke-power of q, we gain 

(mix) = {Cr, Al^_^+ Xm Cm-l}\^^Q, m > 1, < X < +00, (73) 

where 

■^m = Cm - tt 0m,l + ^m,0, (74) 

and the mth-order deformation equation 

VVm(a;, ;2) = 0, m > 1, 2 < 0, < X < +00, (75) 
subject to the boundary condition on the known free surface z = 0: 



r Ij, \ 12 "^Vm , d(j)m 



dx"^ dz 
and the bottom condition 



= Rm{x), < X < +CX), (76) 

z=0 



dz 

where 



0, z = -1, < X < +00, (77) 



Rni{x) = \ C0 A^_i + Xm Sm-1 - Sm\\ , < X < +00, (78) 
L J |2;=0 

Xn = I ?' ^l!'" " - I' (79) 
^ \ 1, when n> 1. ^ ' 

The detailed derivations of A^_^, Sm-i, Sm with all related formulas are given 

explicitly in the Appendix. Note that, unlike Liao & Cheung |15] and Tao et al |48j . 
we explicitly give all formulas in details so that high-order approximations can be 
gained more efficiently. 

Note that the dimensionless phase speed a of the peaked solitary waves is unknown 
up to now. According to the linear wave theory mentioned in § [3l the peaked solitary 
waves exist only when 

2 tan/c , TT 

a = — ; — , niT < k < mv -\ — , (80) 
k 2 

where n > is an integer. In this section, we consider only the primary decay-rate 
< ko < 7t/2, i.e. k = ko. If the above expression also holds for the fully nonlinear 
wave equations, the auxiliary linear operator defined by (1621) has the property 

£{cos[A;(2 + l)]e-'^^} = 0, x>0,k = ko>0, (81) 

and the corresponding inverse operator of C defined by fl76|) has the property 

- 1 , , , cos[nA;(z + 1)1 exp(— nfcx) , , , , 

C-'{exp{-nkx)}= \ \ ^ \ k>0, n^l, x>0, 82 
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where n > 2 is an integer. Note that the above expression does not hold when n = 1. 
Fortunately, it is found that Rm{,x) indeed does not contain the term exp(— fcx) as long 
as the phase speed is given by = tan(/c)//c, where < < 7r/2. Mathematically, 
this is because the nonlinear terms of ( |56ll do not contain the term exp{—kx) at all, 
since 

exp(— m/cx) X exp{—nkx) = e~(™'+")'''' 

with m + n > 2 for any integers m > 1 and n > 1. So do the linear terms of ( 156|1 . 
since 




+ — j ^^&n cos[nk{z + 1)] exp{—nkx) > 

' n=l ) 



^^(nA;) 'ya^{nk) cos{nk) — sin(nA;)] 6„ exp(— nfcx) 

n=l 

+ 00 

^^(nA;) [a'^{nk) cos{nk) — sin(nA;)] 6„exp(— n/cx) 



z=0 



+00 



n=2 



does not contain the term exp{—kx), too. This is the essential reason why the phase 
speed 

a = (83) 

given by the linear wave equations still holds for the fully nonlinear wave equations 
([2]) to ( |TOl) . Physically speaking, the dimensionless phase speed cj \fgD of the peaked 
solitary waves has nothing to do with the wave height: this is quite different from the 
traditional periodic and solitary progressive waves. It indicates that the new kind of 
peaked solitary waves are non- dispersive. We will illustrate this point later. 

Keeping fl83|) in mind and using the property f l82|) of the inverse operator it is 
straightforward to gain the common solution of the high-order deformation equation 

dzg) to dn]): 

(\)^{x,z)=(\r^[x,z)^ A^zo'A[k{\^z)\e-^'', x > 0, (84) 

where 0^(x, z) = L^\Rm{x'^ is a special solution, and the coefficient is deter- 
mined by the given wave height 



m+l 



VlimC„(x) = i/^. (85) 



n=l 



This is mainly because, according to fl75]) . Cm+i(a;) is dependent upon (^^(x, z) that 
contains the unknown parameter A^, where m > 1. Note that, according to fl52]) . 
(prn{x, z) is in the form of (152!) and thus automatically satisfies the Laplace equation ([2]), 
the bottom condition ([5]) and the bounded condition (fTTl) . Thus, using the explicit 
formulas given in the Appendix, it is computationally efficient to gain high-order 
analytic approximations successively, especially by means of computer algebra system 
such as Mathematica and Maple, since our approach needs only algebra computations. 
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For example, using the initial guess ( 15^ and ( !73|) . we directly have 

Ci(x) = 



d(f)o 1 



z=0 



Q;cos(/c)e + ^-e 



X > 0. 



^6) 



Thus, at the first-order of approximation, we have an algebraic equation for the given 
wave height 

Hw = CrjkAo (a cos k + -kAoj , 



which gives two different solutions 



An = k 



-1 



—a cos k ± \/ a'^ cos^(A;) + 2if^/ 



57) 



We simply choice 



Ao = -k-' 



a cos k — J cos^(A;) + 2Hu,/cri 



to calculate Aq for a given if^, since it has a smaller absolute value. 
Furthermore, using the initial guess (1541) . we have 

= kAo {a^k cos k - sin k) e''^^ + 2aA;3^2g-2fcx ^ ^4^3 cos(A;)e-3^^ 

Using the phase speed the term exp(— A;x) of the above expression disappears, 

say. 



= 2ak'Aie-'^'' + cos(A;)e-'''", a; > 0. 
Thus, the first-order deformation equation reads 

V^0i(a;, = 0, 2: < 0, < X < +00, 

subject to the boundary condition on the known free surface z = Q: 



^9) 



£1 



a 



+ 



dx"^ ' dz 
and the bottom condition 

dz 



[2ak^Ale-^''^ + k^Al cos(A;)e-=^'=^] , (90) 



z=0 



2 = —1, < X < -l-cxo. 



(91) 



Using the property of the inverse operator fl82l) . it is easy to gain the common solution 

, ( ak'^Alcos\2k(z+l)]e-^'''' k^Af.coskcos\3k(z + l)]e-^'''' 



2a2A; cos(2A;) - sin(2A;) 3[3a^k cos(3A;) - sin(3fc)] 

+ Aicos[k{z + l)]e''''', 0<x<+oo, (92) 
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where Ai is an unknown constant to be determined. Similarly, using ( !73|) . we gain 
C2{x), which contains the unknown constant Ai. Then, for the given wave height H^, 
we have a linear algebraic equation 

= limCi(a;) + limC2(a;), 

x^O x^O 

which determines Ai. In this way, (j)i{x,z) is completely determined. Similarly, we 
further gain 02 (a;, z), Csl^); cind so on. Finally, using the symmetry (IT^ . we gain the 
wave elevation ({x) and the velocities u{x, z),v{x, z) in the interval —oo < x < and 
< X < +00. At X = 0, the continuous horizontal velocity u(0, z) = U{z) is given by 
iQ, and the vertical velocity v{0,z) = is given directly by the symmetry condition 

Our computations confirm that, for all m > 0, Rm{.x) in (17^ indeed does not 
contain the term exp(— fcx) at all. Thus, the fully nonlinear wave equations ([2]) - (ITT]) 
indeed give the same dimensionless phase speed a = ^Jtax^k/k as that by the linear 
ones. Therefore, the phase speed of the peaked solitary wave has nothing to do with 
the wave height H^, say, the peaked solitary waves are non-dispersive! This is indeed 
completely different from the traditional periodic and solitary waves with smooth 
crest. This unusual characteristic clearly demonstrates the novelty of the new peaked 
solitary surface waves. We will confirm and discuss this interesting characteristic of 
the new kind of peaked solitary waves later. 

Note that our HAM-based analytic approach mentioned above is rather similar 
to those by Liao & Cheung [15] and Tao et al |18] for the traditional progressive 
waves in deep and finite water, except that we use here the symmetry condition f[T^ . 
the evanescent-mode base-function fl2B]) . and besides regard the dimensionless phase 
speed a as a given constant. 

Finally, we should emphasize that, unlike perturbation methods, our HAM-based 
analytic approach does not need any assumptions about small/large physical parame- 
ters. More importantly, both of (f){x,z) and C,{x) contain the two convergence-control 
parameters and c^, which provide us a convenient way to guarantee the convergence 
of approximation series, as illustrated below. 

4.2 Convergence of series solution 

Note that, unlike perturbation results, 4>m{x, z) and Cm{,x) gained in above-mentioned 
analytic approach contain two convergence-control parameters and Cj^, which pro- 
vide us a convenient way to guarantee the convergence of the series (1711 ) and (1721) , as 
shown below. Obviously, the convergence rate of the series fITT] ) and fl72l) is greatly 
influenced by and Cr,. As pointed out by Liao & Cheung [45] and Tao et al |18] . 
one can choose c<^ = — 1 and c,, = — 1 for weakly nonlinear waves. 

First, let us consider the case of /c = 1 and Hu, = 1/20, with the corresponding 
dimensionless phase velocity a = c/y/gD = A/tan(l) = 1.24796. Since the wave 
height is only 5% of the water depth D, the nonlinearity is weak. Thus, as suggested 
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Order of approx. 


U{-1) 


f/(-0.5) 


[/(-0.25) 




C'(o+) 


1 


0.07222 


0.06570 


0.05762 


0.04289 


-0.04690 


3 


0.06833 


0.06236 


0.05466 


0.04205 


-0.04859 


5 


0.06796 


0.06219 


0.05489 


0.04213 


-0.04823 


10 


0.06799 


0.06221 


0.05490 


0.04213 


-0.04823 


15 


0.06799 


0.06221 


0.05490 


0.04213 


-0.04823 


20 


0.06799 


0.06221 


0.05490 


0.04213 


-0.04823 


25 


0.06799 


0.06221 


0.05490 


0.04213 


-0.04823 



Table 1: Analytic approximations of U{z) — u{0,z) and C'(0+) in the case of Hyj — 
1/20 and A; = 1 by means of — —1 and = —1. 



Order of approx. 


U{-1) 


f/(-0.5) 


f/(-0.25) 




C'(o+) 


1 


-0.07047 


-0.06377 


-0.05101 


-0.03840 


0.05248 


3 


-0.08133 


-0.06737 


-0.05218 


-0.03779 


0.05220 


5 


-0.08140 


-0.06749 


-0.05218 


-0.03772 


0.05192 


10 


-0.08145 


-0.06750 


-0.05218 


-0.03772 


0.05183 


20 


-0.08145 


-0.06750 


-0.05218 


-0.03772 


0.05183 


25 


-0.08145 


-0.06750 


-0.05218 


-0.03772 


0.05183 



Table 2: Analytic approximations of U{z) — u{0,z) and C'(O-i-) the case of Hyj — 
—1/20 and A; = 1 by means of = —1 and = —1. 
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Order of approx. 


U(-l) 


f/(-0.5) 


f/(-0.25) 






1 


0.1696 


0.1543 


0.1347 


0.09362 


-0.08561 


3 


0.1246 


0.1180 


0.1090 


0.08837 


-0.09332 


5 


0.1270 


0.1196 


0.1098 


0.08813 


-0.09142 


10 


0.1254 


0.1183 


0.1090 


0.08788 


-0.09285 


15 


0.1254 


0.1183 


0.1090 


0.08789 


-0.09299 


20 


0.1254 


0.1183 


0.1090 


0.08789 


-0.09299 


25 


0.1254 


0.1183 


0.1090 


0.08789 


-0.09299 



Table 3: Analytic approximations of U{z) = u{0,z) and C'(0+) the case of = 
1/10 and A; = 1 by means of = —1/2 and Cr^ = —1. 



by Liao & Cheung [15] and Tao et al [H], we choose = — 1 and c,, = — 1 for 
such a kind of weakly nonlinear wave problem. It is found that, the corresponding 
series of analytic approximation indeed converges quickly, as shown for examples in 
Table [T]for C'(0+) the horizontal velocity U{z) = u{0, z) beneath the wave crest at 
z = —1, z = —1/2, z = —1/4 and z = H^, respectively, where 0+ denotes x — )■ from 
the right along the x axis. It is found that the velocity potential z) converges 
quickly in the whole domain x G (0, +oo) and z < ({x), as shown for examples in 
Fig. [3] for the corresponding horizontal velocity profile U{z) = u{0,z) beneath the 
crest. This confirms that the new kind of peaked solitary waves is indeed a solution 
of the fully nonlinear wave equations ([2]) - ffTTl) ! 

Secondly, let us consider the case with k = 1 and H^j = —1/20, with the same 
dimensionless phase velocity a = c/y^gD = 1.24796. It is found that the correspond- 
ing series of analytic approximations given by = — 1 and = — 1 converge quickly 
in the whole domain x > 0, as shown for examples in Table [2] for C'(0+) ^ind the hor- 
izontal velocity U{z) = u{0, z) beneath the crest aX z = —1, —0.5, —0.25 and z = H^, 
respectively. Besides, the corresponding velocity potential z) converges quickly 
in the whole domain x G (0, +oo) and z < C{x), as shown for examples in Fig. IHfor 
the horizontal velocity profile U{z) = m(0, z) beneath the crest. This confirms that 
the new kind of peaked solitary waves in the form of depression is truly a solution of 
the fully nonlinear wave equations ([2]) - (fTTj) ! 

Similarly, in the case of = 1/2 and H^j = ±1/20, with the corresponding di- 
mensionless phase speed a = c/ yJgD = 1.04528, we gain convergent series of analytic 
approximations of (f>{x,z) and ({x) in the whole domain x E (0, +oo) and z < ({x) 
by means of = — 1 and = —1, respectively. 

Furthermore, let us consider the case of = 1 and Hyj = 0.1, with the corre- 
sponding dimensionless phase velocity a = cj^fgD = 1.24796. Since the wave weight 
increases to 10% of the water depth, the nonlinearity becomes stronger. As suggested 
by Liao & Cheung |15] and Tao et al [H], we should choose convergence-control pa- 
rameters and with smaller absolute values for the higher nonlinearity. It is found 
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Order of approx. 


f/(-l) 


f/(-0.5) 


t/(-0.25) 




C'(o+) 


1 


-0.1418 


-0.1191 


-0.09328 


-0.07474 


0.1091 


3 


-0.1704 


-0.1343 


-0.09684 


-0.07211 


0.1130 


5 


-0.1768 


-0.1368 


-0.09664 


-0.07070 


0.1107 


10 


-0.1801 


-0.1379 


-0.09650 


-0.07016 


0.1079 


15 


-0.1805 


-0.1380 


-0.09649 


-0.07013 


0.1075 


20 


-0.1806 


-0.1380 


-0.09648 


-0.07012 


0.1075 


25 


-0.1806 


-0.1380 


-0.09648 


-0.07012 


0.1075 



Table 4: Analytic approximations of U{z) = u{0,z) and C'(0+) the case of = 
—1/10 and A; = 1 by means of C(j, = —3/4 and c,, = —1. 



0.1 r 




0.03 0.04 0.05 0.06 0.07 0.08 

U(z) 

Figure 3: Analytic approximations of the dimensionless horizontal velocity profile 
U{z) = m(0, z) beneath the crest in the case of A; = 1 and = 0.05 given by = —1 
and = —1. Dashed-line: zeroth-order of approx.; Dash-dotted line: Ist-order of 
approx.; Solid line: 4th-order of approx.; Symbols: 25th-order of approximation. 
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U(z) 

Figure 4: Analytic approximations of the dimensionless horizontal velocity profile 
U{z) = u{0, z) beneath the crest in the case of A; = 1 and Hyj = —0.05 given by = 
— 1 and Crj = —1. Dashed-line: zeroth-order of approx.; Dash-dotted line: Ist-order of 
approx.; Solid line: 4th-order of approx.; Symbols: 25th-order of approximation. 




U(z) 

Figure 5: Analytic approximations of the dimensionless horizontal velocity profile 
U{z) = m(0, z) beneath the crest in the case of = 1 and = 0.1 given by c<^ = —1/2 
and = —1. Dashed-line: zeroth-order of approx.; Dash-dotted line: 2nd-order of 
approx.; Solid line: 6th-order of approx.; Symbols: 25th-order of approximation. 
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Figure 6: Analytic approximations of the dimensionless horizontal velocity profile 
U{z) = m(0, z) beneath the crest in the case of = 1 and H^j = —0.1 given by = 
—3/4 and = —1. Dashed-line: zeroth-order of approx.; Dash-dotted line: Ist-order 
of approx.; Solid line: lOth-order of approx.; Symbols: 25th-order of approximation. 




Figure 7: Analytic approximations of elevation of the peaked solitary waves when 
k = 1 (corresponding to c/ yJgD = 1.24796). Solid line: 5th-order approximation 
when Hy^ = 0.1 given by c,^ = —0.5 and = —1; Filled circles: 25th-order approx- 
imation when Hyi; = 0.1 given by = —0.5 and = —1; Dashed line: 5th-order 
approximation when Hy^ = 0.05 given by = —1 and = —1; Open circles: 25th- 
order approximation when = 0.05 given by = — 1 and = —1. 
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X 



-QI 




Figure 8: Analytic approximations of elevation of the peaked solitary waves when 
k = 1 (corresponding to c/y/gD = 1.24796). Solid line: 5th-order approximation 
when = —0.1 given by = —0.75 and = —1; Filled circles: 25th-order 
approximation when = —0.1 given by = —0.75 and = —1; Dashed hne: 
5th-order approximation when = —0.05 given by = —1 and = —1; Open 
circles: 25th-order approximation when Hy^ = —0.05 given by = —1 and = —1. 




Figure 9: Analytic approximations of ({x) of the peaked solitary waves when = 
0.05 by means of = —1 and = —1. Solid line: 5th-order approximation when k = 
1/2 (corresponding to cj^fgD = 1.04528); Filled circles: 25th-order approximation 
when k = 1/2; Dashed line: 5th-order approximation when k = 1 (corresponding to 
c/y/gD = 1.24796); Open circles: 25th-order approximation when k = 1. 
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Figure 10: Analytic approximations of ({x) of the peaked solitary waves when 
Hw = —0.05 by means of = —1 and = —1. Solid line: 5th-order approxi- 
mation when k = 1/2 (corresponding to cj^fgD = 1.04528); Filled circles: 25th-order 
approximation when k = 1/2; Dashed line: 5th-order approximation when k = 1 
(corresponding to cj \fgD = 1.24796; Open circles: 25th-order approximation when 
k = l. 

that the series of analytic approximations given by = —1/2 and = — 1 converges 
quickly, as shown in Tableland Fig.[5]for the horizontal velocity profile U{z) = u{0, z) 
beneath the wave crest. Similarly, in the case of = 1 and Hyj = —0.1, we gain con- 
vergent series of analytic approximation by means of = —3/4 and = —1, as 
shown in Table H] and Fig. [HI This illustrates that the two convergence- control pa- 
rameters and Crj indeed provide us a convenient way to guarantee the convergence 
of approximation series. Note that the absolute value of the horizontal velocity at 
bottom in the case of Hy^ = —0.1 is 44% larger than that in the case of Hy^ = 0.1. So, 
there does not exist symmetry between these two wave elevations for Hy^ = 0.1 and 
Hw = —0.1, respectively. 

It should be emphasized that, in the case of /c = 1, we gain convergent series 
solutions of the peaked solitary waves with the same phase speed but different positive 
and negative values of Hyj such as Hyj = ±0.05 and = ±0.1. This confirms that the 
phase speed of the peaked solitary waves indeed does not depend upon the wave height, 
say, the peaked solitary waves are non- dispersive. This is an unusual characteristic of 
the peaked solitary waves. 

Finally, using the symmetry f lT^ . it is straightforward to gain the wave elevation 
in the whole interval — oo < x < ±oo. As shown in Figs. [7]to[T0l the wave elevation 
({x) also converges quickly in all of above-mentioned cases. In case of k = 1, the wave 
elevations when Hyj = ±0.1 are compared with those with Hyj = ±0.05, as shown in 
Figs. [7]and[Hl It is found that, for the same value of k, the larger the value of \Hy,\, 
the faster ({x) decays to zero. Note also that the wave elevation ({x) with larger k 
decays to more quickly, as shown in Figs. |9] and [101 In other words, the larger the 
value of k, the faster C(x) — )■ 0. This provides us a physical meaning of the parameter 
k. For this reason, we call k the decaying-rate parameter. 

Note that our HAM-based analytic approach contains two convergence-control 
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parameters c^f, and c,,. As illustrated above, we can gain convergent series of analytic 
approximations by choosing proper values of and c^, which indeed provide us a 
convenient way to guarantee the convergence of approximation series. 

Note that, as proved in general by Liao [33l|36], each series solution given by 
the HAM satisfies its original equations as long as it is convergent. So, all of these 
convergent series of z) and C,{x) are solutions of the fully nonlinear wave equations 
(E]) - ( ITTj) . as further confirmed below. 



4.3 Validation check of analytic approximations 



Note that the velocity potential z) is expressed in the form flS^ . which automat- 
ically satisfies the Laplace equation (|2]), the bottom condition ([5]), and the bounded 
condition (fTT]) . Thus, it is only necessary for us to check the two nonlinear boundary 
conditions ([3]) and (j4]) which are satisfied on the unknown wave elevation C,{^)- 

To check the validation of our analytic approximations, we define the averaged 
residual squares of the two free surface boundary conditions 



1 ^ 



n=l 
M 



-E 



n=l 



(93) 
(94) 



where 



n=0 



n=l 



are the mth-order approximation of (j){x,z) and C{x), respectively, and 



n 



l<n<M, 



with large enough xr and M. For all results given below, we choose x^ = 10 and 
M = 100, if not mentioned. Since the potential velocity 0(x, z) and the wave elevation 
({x) decay exponentially, xr = 10 is large enough. 

In case of = 1 and Hw = ±0.05, the averaged residual squares £^ and £^ of 
the corresponding analytic approximations obtained by c<^ = — 1 and = — 1 decay 
quickly to the level 10~^^ as the order of approximation increases to 25, as shown in 
Table O In other words, our 25th-order approximation of (j){x,z) and ({x) satisfies 
the Laplace equation ([2]), the bottom condition and the bounded condition flTTl) 
exactly, and besides the two nonlinear free surface boundary conditions ([3]) and (jlj) 
very accurately (to the level 10~^^). In addition, u{0,z) determined by ([9]) converges 
quickly so that U (z) is uniquely determined. Therefore, without doubt, our convergent 
analytic approximation is a very accurate solution of the generalized, fully nonlinear 
wave equations ([2]) to (ITT]) . Similarly, £^ and £^ decays to the level 10~^^ in the 
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Order of approx. 


Hw = 


0.05 




Hw = 


-0.05 




m 








f<t> 








1 


6.59 xlQ-^ 


9.21 xlO 


-9 


2.25 xlO- 


-6 


9.39 xlO 


-9 


3 


1.40 xlO-s 


1.57 xlO- 


-9 


8.09 xlO- 


-9 


5.90 xlO- 


-10 


5 


3.32 xlO-ii 


2.90 xlO- 


-12 


7.80 xlO- 


11 


1.53 xlO- 


-11 


10 


9.71 xlO-15 


7.00 xlO' 


-16 


8.96 xlO- 


16 


2.42 xlO- 


-16 


15 


7.33 X 10^19 


1.68 xlO- 


-19 


6.97 xlO" 


20 


1.05 xlO- 


-20 


20 


4.43 X 10^22 


4.83 xlO- 


-23 


2.40 xlO" 


24 


6.45 xlO" 


-25 


25 


2.23 xl0~25 


2.09 xlO- 


-27 


4.56 xlO- 


28 


3.96 xlO- 


-29 



Table 5: Averaged residual squares of the two nonlinear free boundary conditions ([3]) 
and (jl]) in the case of k = 1 and H^ = ±0.05 by means of = —1 and = —1, with 
the corresponding dimensionless phase speed cj \fgD = 1.24796. 



case of = 1 and Hw = ±0.1, and to the level 10~^^ in the case of k = 1/2 and 
Hw = ±0.05, respectively, as shown in Tables M and [71 Theses guarantee that the 
corresponding analytic approximations of (f){x,z) and ({x) are indeed quite accurate 
solutions of the fully nonlinear wave equations ([2]) to ( ITTj) . respectively. All of these 
confirm once again the mathematical proof of Liao [33ll36] that each convergent series 
solution given by the HAM satisfies its original equations in general. 

It should be emphasized that, in the case of k = 1, correspond to the same 
dimensionless phase speed c/y/gD = 1.24796, we gain very accurate solutions of 
solitary waves with different positive and negative wave heights such as Hw = ±0.05 
and Hw = ±0.1, respectively. This confirms that the phase speed of the peaked 
solitary waves is indeed independent of the wave height Hw- 

In fact, one can choose the optimal values of c,^ and by the minimum of 

^tic^^c^) and £^{c^,Cr,), say. 

It is found that, by means of the optimal values of c,^ and c^, the corresponding series 
of analytic approximations often converge more quickly. 

In addition, we also apply perturbation method to check the validity of our ana- 
lytic solutions. It is found that the first-order perturbation approximation of z) 
is exact the same as f l92p given by our HAM-based approach in the case of = — 1 
and = —1. This confirms the validity of our analytic approximations in a different 
way. 

All of these demonstrate that the convergent series obtained by our HAM-based 
approach are indeed the solutions of the generalized, fully nonlinear wave equations 

to dm). 
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Order of approx. 




= 0.1 




Hw = 


-0.1 




m 


(c^ = -0.5, = -1) 




(c0 = -0.75, = -1) 






c<t> 






f<l> 






1 


1.48 xlO"* 


7.21 xlO- 


-7 


5.89 xlO"^ 


3.67 XlO- 


-7 


3 


1.63 xlO~^ 


8.63 xlO 


-8 


5.84 xlO-^ 


1.21 xlO- 


-7 


5 


1.39 xlO-^ 


3.96 XlO- 


10 


6.83 xlO-^ 


1.11 xlO- 


-8 


10 


5.96 xlO-i° 


3.09 xlO- 


-11 


8.63 xlO-9 


2.31 xlO- 


10 


15 


1.30 xlO-12 


3.70 xlO- 


-14 


3.29 xlO'i^ 


1.95 xlO- 


12 


20 


4.34 X 10-^3 


1.88 xlO" 


-15 


1.17 xlO-^2 


3.69 xlO" 


14 


25 


2.25 xlO-13 


4.52 xlO" 


-16 


2.51 xlO-i'^ 


8.11 xlO- 


16 



Table 6: Averaged residual squares of the two nonlinear free boundary conditions ([3l) 
and (jl]) in the case of = 1 and = ±0.1, with the corresponding dimensionless 
phase speed c/ ^/gD = 1.24796. 



Order of approx. 


Hw = 


0.05 




Hw = 


-0.05 




m 


C(f> 
'-'m 






p<f> 








1 


2.75 xlO"^ 


1.04 xlO- 


-6 


3.74 XlO" 


-7 


9.09 XlO 


-7 


3 


3.76 xlO-i° 


4.85 xlO 


-8 


4.48 xlO- 


-10 


9.61 xlO 


-9 


5 


4.13 xlO-12 


2.00 xlO- 


-9 


4.85 xlO- 


-13 


4.18 xlO- 


-12 


10 


1.80 xlO-i^ 


2.21 xlO- 


-12 


1.08 xlO- 


-16 


3.16 xlO- 


-17 


15 


5.23 X 10^15 


1.55 xlO- 


-15 


7.93 xlO- 


-20 


1.06 xlO- 


-19 


20 


1.86 xlO-i^ 


1.50 xlO- 


-16 


1.67 xlO- 


-23 


1.52 xlO' 


-24 


25 


2.56 xlO-is 


9.46 xlO" 


-18 


3.87 XlO- 


-29 


2.57 xlO- 


-28 



Table 7: Averaged residual squares of the two nonlinear free boundary conditions (El) 
and (jl]) in the case of /c = 1/2 and H^ = ±0.05 by means of = —1 and Cr, = —I, 
with the corresponding dimensionless phase speed c/y/gD = 1.04528. 



34 



4.4 Characteristics of peaked solitary surface waves 

Based on the so-called evanescent-mode base-functions fl26l) . the new kind of solitary 
waves have some unusual characteristics that are quite different from those of the 
traditional waves with smooth crest. 

First, the new kind of solitary waves have a peaked wave crest, since C'{x) is 
discontinuous at x = 0, i.e. C'(0+) 7^ C'(O-), where 0+ and 0_ denote x — )■ from 
the right and left along the x axis, respectively. For example, in the case of = 0.1 
and k = l, ('(0+) = -0.09299, but C'(O-) = 0.09299. This is obviously different from 
traditional periodic and solitary waves which are infinitely differentiable. 

Secondly, the new kind of peaked solitary waves may be in the form of depression, 
which has been reported for internal waves but never for surface ones. Mathematically, 
it is straightforward to gain such kind of solitary waves of depression even by means 
of the linear wave equations, as shown in § [31 

Third, unlike traditional periodic and solitary waves, the dimensionless phase 
speed of the new kind of peaked solitary waves depends only upon k, the so-called 
decaying-rate parameter, but has nothing to do with the wave height H^. So, in the 
same water depth D, the new kind of solitary waves with the same k but different wave 
height Hyj may propagate with the same phase speed, where if^ may be either positive 
or negative. For example, it is found that, in the case of /c = 1, all of the peaked 
solitary waves with Huj = ±0.1 or H^j = ±0.5 propagate with the same phase speed 
c = 1.24796a/^^ : in these cases, we gain different convergent series solution with the 
same phase speed, as shown in § 14.21 On the other side, the peaked solitary waves 
with the same wave height but different decay-rate parameter k may propagate 
with different phase speed! These are quite different from the traditional periodic and 
solitary waves whose phase speed strongly depends upon wave height. In other words, 
the traditional periodic and solitary waves with smooth crest are dispersive, but the 
new peaked solitary waves are non- dispersive. 

Finally, as shown in Tables [T] to H] and Figs. 2 to |6l the horizontal bottom 
velocity beneath the crest of the new kind of peaked solitary waves is always larger 
than that at crest. For example, in the case of A; = 1 and Hu, = 0.1, the horizontal 
velocity at bottom beneath the crest is 43% larger than that at crest, as shown in 
Table [31 Especially, as shown in Table [H in the case of k = 1 and = —0.1, the 
horizontal velocity at bottom beneath the crest is even 158% larger than that at crest! 
In general, for the same x, the horizontal velocity u{x, —1) (at bottom) has always 
a larger absolute value than u{x,({x)) on the free surface, as shown for example in 
Figs. [TTland[T2l This is quite different from the traditional periodic and solitary waves 
whose horizontal velocity at bottom is always less than that on surface. Besides, these 
also illustrate that there does not exist symmetry between any two kinds of peaked 
solitary waves with Hy^ > and H^^ < 0, respectively. 

Furthermore, as shown in Figs. [7] and [H in the case of the same k, the peaked 
solitary wave with larger value of \Hw\ is sharper at crest. It is found that, in general, 
the larger the value of \Hyj\, the sharper the peaked solitary waves at crest, as shown 
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Figure 11: Horizontal velocity at bottom and on free surface when k = 1 and 
Huj = 0.05 by means of c<^ = —1 and = —1. Solid line: 3rd-order approx. of 
u{x,—l) (at bottom); Dashed line: 3rd-order approx. of u{x,({x)) (on free surface); 
Symbols: the corresponding 25th-order approximations. 




Figure 12: Horizontal velocity at bottom and on free surface when k = 1 and 
Hw = —0.05 by means of = —1 and = —1. Solid line: 3rd-order approx. of 
m(x, — 1) (at bottom); Dashed line: 3rd-order approx. of u{x,({x)) (on free surface); 
Symbols: the corresponding 25th-order approximations. 
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Figure 14: 
line: if,,, = 



C'(0+) (as X — )■ from right) versus k. Solid line: H^j = 0.05; Dashed 
-0.05. 
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in Fig. [131 Note that all of these peaked solitary waves with the same k but different 
Hyj propagate with the same phase speed! In addition, as shown in Figs. [9] and [TOl 
in the case of the same value of if^, the larger the value of fc, the sharper the peaked 
solitary wave at crest. This also holds in general: the larger the value of the sharper 
the peaked solitary waves at crest, as shown in FigHH Therefore, generally speaking, 
the larger the values of k and the sharper the peaked solitary waves. 

All of these unusual characteristics clearly indicate the novelty of the peaked 
solitary waves. It should be emphasized that these so-called peaked solitary waves 
given by the linear wave equations in § [3] also have the same unusual characteristics. 

Finally, it should be emphasized that all results given in § H] are gained in the 
case of the primary decay-rate k = k^, where < fco < vr/2 is determined by the 
transcendental equation (1281) . In this case, the velocity potential (f){x,z) and wave 
elevation ({x) governed by the generalized, fully nonlinear wave equations ([2]) to (ITT!) 
can be indeed expressed by (!52|) and (1531) . respectively, as illustrated in §|11 However, 
as mentioned in § [31 for a given a, the transcendental equation fl28p has an infinite 
number of solutions k^,, where z/vr < k^ < z/vr + 7r/2 for an integer u > 0. So, in 
a similar way, we could gain convergent series approximations of the fully nonlinear 
wave equations for other decay-rate k^, and should obtain the same conclusions qual- 
itatively. Obviously, it becomes more complicated if the wave elevation is expressed 
by evanescent-mode with more than two different decay-rates. For example, in the 
case of two different decay-rates /cq and ki, the wave elevation should be in the form 

+ 00 +00 

({x) = g^.o exp(— mfco|x|) + y~^ao^„exp(— nfci|a:|) 

m=l n=l 

+ 00 +00 

+ '^'^am,neyip[-{mko + nki)\x\], (96) 

m=l n=l 

which is much more complicated. Considering the length of this manuscript, we will 
show these more complicated peaked solitary waves in details in the near future. 

5 Concluding remarks, discussions and some the- 
oretical predictions 

In this article, the progressive surface waves with permanent form are described in 
a generalized way by means of the PDEs ([2]) to flTTl) in the interval < x < +oo, 
whose solution is then extended to the interval (— oo, 0) by means of the symmetry 
f[T2|) . Besides, at a; = 0, the vertical velocity v{0,z) = is directly determined by 
the symmetry f[T^ . the horizontal velocity u{0,z) is determined by dHj) as a limit of 
u{x, z) when x — )■ 0, and the wave elevation ^(0) is determined by f[TUl) as a limit of 
({x) when x — )■ 0, respectively. As shown in § 2, this approach admits the traditional 
progressive waves with smooth crest. However, as shown in § 3 and § 4, it also admits 
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the peaked solitary waves, which have many unusual characteristics. Therefore, this 
approach is more general than the traditional ones. 

By means of the so-called evanescent-mode base-functions fl26l) . a new type of 
peaked solitary surface waves is found first by using the linear wave equations in § [3] 
and then confirmed by the generalized, fully nonlinear wave equations (E]) to ( |TT|) in 
§ m Using the homotopy analysis method (HAM), we gain the convergent series of 
the velocity potential (j){x, z) and the wave elevation ({x) for different primary decay- 
rate parameter k and positive/negative wave height by means of the evanescent- 
mode base-functions (!26l) . The validity of these convergent analytic approximations 
are carefully checked: all linear governing equations and boundary conditions are 
automatically satisfied, and besides the two nonlinear free surface boundary conditions 
are satisfied very accurately, as shown in Tables [5] to [3 So, we are quite sure that our 
convergent series approximations of the peaked solitary waves given in § H] are indeed 
the solutions of the generalized, fully nonlinear wave equations ([2]) to ffTTj) . 

It is found that these peaked solitary waves expressed by the evanescent-mode 
base-functions fl26p have many unusual characteristics different from the traditional 
periodic and solitary ones expressed by the smooth base-functions f|T6l) . First, it has 
a peaked crest. Secondly, it may be in the form of depression, corresponding to a 
negative wave height H^, which has been often reported for interfacial solitary waves 
but never for free-surface solitary ones, to the best of author's knowledge. Third, its 
phase speed has nothing to do with wave weight H^, say, the peaked solitary waves 
are non-dispersive. Finally, its horizontal velocity at bottom is always larger than 
that on free surface. All of these are so different from the traditional periodic and 
solitary waves with smooth crest that they clearly indicate the novelty of the peaked 
solitary ones reported in this article. 

All of these unusual characteristics come from the so-called evanescent-mode base- 
functions f l2B]) of the peaked solitary waves, which are essentially different from the 
smooth base-functions (fT6|) of the traditional periodic and solitary waves, although 
both of them automatically satisfy the Laplace equation ([2]), the bottom condition 
([5]) and the bounded condition (fTTI) . It should be emphasized that both of them 
are widely used and can be found in the textbook [2] of water waves, although the 
smooth base-functions (1161) are currently more familiar and thus might be regarded as 
the mainstream. Note that the traditional base-functions (I16p are infinitely differen- 
tiable everywhere, but the derivatives of the evanescent-mode base-functions fl26|) with 
respect to x are not differentiable at a; = 0. This essential difference of the two kinds 
of base-functions f|T6|) and f l26|) is the origin of the completely different characteristics 
between the smooth and peaked waves. 

Note that all traditional smooth periodic and/or solitary progressive waves can 
be derived by means of the generalized, fully nonlinear wave equations (|2]) to f|TT]) 
with the symmetry f lT2|) . In other words, Eqs. ([2]) to f lTTj) with the symmetry f|T2|) 
is consistent with the traditional wave theory. Note that the peaked solitary waves 
are lost if the Laplace equation ([2]) is defined in the whole interval (— oo, +oo). Thus, 
the generalized, fully nonlinear wave equations ([2]) to f ITT]) with the symmetry f[T2]) 
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provide us a more general approach to investigate the progressive waves. 

It is found that the generahzed, fully nonlinear wave equations ([2]) to ffTTl) contain 
two different types of solutions: one (the traditional progressive periodic and solitary 
wave) has an infinitely differentiable wave elevation ({x) with ^'(0) = 0, the other (the 
peaked solitary wave reported in this paper) has a peaked crest with the discontinuous 
C'(0), corresponding to the two different base-functions (fT6|l and ( !26|) . respectively. It 
should be emphasized that, according to the theory of differential equations, the 
Laplace equation (|2]) needs only one boundary condition ([9]) at x = 0, so that any 
other smoothness conditions at x = such as the infinitely differentiable velocity 
potential is unnecessary and thus should be avoided, since they may lead to the loss 
of the peaked solitary waves. 

In theory, such kind of gravity surface waves with peaked crest are not new at all. 
It is well-known that the limiting gravity wave has a corner crest with 120 degree, as 
pointed out by Stokes [1] in 1894. In 1993, Camassa & Holm [5] found the peaked 
solitary waves by means of the CH equation ([T]) in the special case of u = 0, which is 
the same as (I35p found in this article by means of the linear wave equations. Currently, 
it is found by Liao [11] that the CH equation ([1]) also admits peaked solitary waves 
even in the case of cu 7^ 0. Besides, the closed-form solutions of the peaked solitary 
waves of the KdV equation |2] , the Boussinesq equation [1], the BBM equation [3] and 
modified KdV equation are currently found by Liao [10]. Therefore, in theory, nearly 
all mainstream models of shallow water waves admit the peaked solitary waves. Note 
that all of these equations are approximations of the fully nonlinear wave equations 
in shallow water, so that our peaked solitary waves derived from the fully nonlinear 
wave equations ([2]) to (fT2!) well explain why all of these shallow water equations admit 
peaked solitary waves. As mentioned in § 1, it was an open question whether or not 
the fully nonlinear wave equations admit the peaked solitary waves. In this article, a 
positive answer to this open equation is given. 

In practice, it is well-known that solutions related to dam break and shock waves 
are discontinuous. Such kind of discontinuous problems belong to the so-called Rie- 
mann problems [521455] . a classic field of fluid mechanics. Such kind of discontinuity 
(or singularity) of solutions of water wave equations have clear physical meanings, 
which are often solved in different sub-domains by numerical and analytic methods. 
For such kind of Riemann problems, it is unnecessary for the considered differential 
equations and boundary conditions to be satisfied at points with discontinuity. The 
numerical and analytic results of many Riemann problems agree well with experi- 
mental results, indicating that such kind of discontinuity (or singularity) should be 
reasonable not only in mathematics but also in physics. So, physically, the peaked 
solitary waves are not strange at all even from the traditional view-points. Therefore, 
the new type of peaked solitary waves should be reasonable and acceptable, too. 

Due to the symmetry f|T2|) . the vertical velocity at a; = always equals to zero. 
The vertical velocity of the traditional progressive waves with smooth crest is always 
zero at X = 0, and thus is always continuous at crest. However, the vertical velocity 
of the peaked solitary waves expressed by the evanescent-mode base-functions f l2B|) is 
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Figure 15: Possible complicated elevation of three peaked solitary waves with the 
same phase speed and the permanent form ({x) = e~'^'/10 — e~'^~^'/20 + 3e~'^'^^'/100, 
predicted by the linear theory in § [31 

discontinuous at x = 0: it changes sign when we cross the interface x = 0. It should be 
emphasized that such kind of discontinuity of velocity exists even for the traditional 
periodic progressive waves expressed by the smooth base-functions f lT6|) : "the tangen- 
tial velocity changes sign as we cross the surface" of interfacial waves, and "in reality 
the discontinuity, if it could ever be originated, would be immediately abolished by 
viscosity", as mentioned by Lamb [30J (§231, page 371). It should be emphasized 
that both of the smooth base-functions (fT6|) and the evanescent-mode base-functions 
( 126|) are widely used in hydrodynamics and can be found in the textbook [H] of wave 
theories. Therefore, theoretically speaking, they are consistent with each other. So, 
the discontinuity of the vertical velocity of the peaked solitary waves beneath the crest 
is acceptable even from the viewpoint of traditional wave theory. 

On the other hand, such kind of discontinuity (or singularity) of the peaked 
solitary waves can be removed in the following way. Let U (z) denote the horizontal 
velocity at x = 0, corresponding to a progressive wave governed by the fully nonlinear 
wave equations ([2]) to ( ITTj) . Assume that, in the frame moving with the solitary wave, 
one can instantaneously replace the boundary x = by a porous vertical plate, and at 
the same time enforce a horizontal velocity U{z) = u{0, z) through the porous plate. 
Then, the corresponding velocity potential and wave elevation C^{x) in the domain 
X G (0, +oo) are governed by the same wave equations ([2]) to ffTTl) . Therefore, for 
a properly given U{z), one gains either the traditional periodic and solitary waves 
with smooth crest (in the domain < x < +oo), if the smooth base-functions f|T6l) 
are used, or the new type of peaked solitary waves reported in this article, if the 
evanescent-mode base-functions fl2B]) are used, respectively. In this case, the Laplace 
equation ([2]) is satisfied in the whole domain x G (0, +oo) and z < ({x) so that no 
discontinuity (or singularity) exists at all. If necessary, using the symmetry (fT2l) . it is 
straightforward to gain the wave elevation ({x) in the whole domain — oo < x < +oo. 

The peaked solitary surface waves may provide us new explanations for a few nat- 
ural phenomenon and some theoretical predictions. For example, the peaked solitary 
wave has an unusual and interesting characteristic: its phase speed has nothing to do 
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with the wave height Hyj that may be either positive or negative. According to the 
transcendental equation fl28l) . the peaked sohtary waves with small wave height 
may propagate very quickly, since tan(A;)/A; — )■ +00 as k — )■ nvr + 7r/2 for an integer 
n > 0. Thus, all of these peaked solitary waves with small wave height and different 
phase speed may create a huge solitary surface wave somewhere: this gives a new theo- 
retical explanation about the so-called rogue wave that can suddenly appear on ocean 
even when "the weather was good, with clear skies and glassy swells", as reported 
by Graham [56] and mentioned by Kharif [57]. On the other hand, several peaked 
solitary waves with different wave heights if^ and different positions of crest but the 
same decay-rate k can propagate together with a permanent form and the same phase 
speed, as shown for example in Fig. [151 However, such kind of complicated solitary 
waves have never been observed in experiment and practice: this provides us a new 
theoretical prediction. 

In addition, the new type of peaked solitary waves may be in the form of de- 
pression, corresponding to a negative wave height if^. To the best of the author's 
knowledge, such kind of solitary waves of depression have been reported only for in- 
terfacial waves but never for surface gravity waves. Obviously, such kind of peaked 
solitary surface waves of depression should be more difficult to create than the tra- 
ditional ones. However, if this theoretical prediction is indeed physically correct, we 
should observe it in laboratory experiments and/or in practice, sooner or later. This 
is an interesting but challenging work: it could enrich and deepen our understanding 
about solitary waves, no matter the final conclusions are positive or not. 

Possibly, the new kind of peaked solitary surface waves might change some of our 
traditional view-points. For example, solitary waves are often regarded as a nonlinear 
phenomenon. However, we illustrate here, for the first time, that solitary surface waves 
may exist even in a system of linear differential equations! Besides, it is also widely 
believed that solitary waves exist only in shallow water. But, we indicates in this 
article that solitary waves can exist even in a finite water depth, say, D is unnecessary 
to be small. For instance, in the case oi D = 100 meter, k = 1 and the dimensionless 
wave height = 0.05, the corresponding new type of peaked solitary surface wave 
propagates with the 5 meter wave-height in the phase speed c = 1.24796a/^I7 ^ 39.1 
meter per second, which is not very dangerous. However, in the case of D = 1000 
meter, k = 1 and if^, = 0.05, the corresponding peaked sohtary wave propagates with 
the 50 meter wave-height and the phase speed c = 123.5 meter per second, which is 
destructive if it indeed could occur on the earth. Here, it may be worth mentioning 
that, the peaked solitary wave may contain large kinematic energy only in a small 
domain near a; = 0, since the largest horizontal velocity is near the bottom and besides 
the velocity decays exponentially in the x direction. So, there exists mass movement 
only in a small domain near x = 0. Certainly, the mass movement of this peaked 
solitary waves is larger than that of the traditional solitary waves with smooth crest, 
but it might be realistic. Thus, the peaked solitary waves are more dangerous than 
the traditional solitary ones with smooth crest. 

According to the traditional wave theories with smooth crest, the velocity of fluid 
decreases exponentially in the vertical direction (from surface to bottom) so that a 
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submarine far enough beneath ocean surface is safe even if there are huge waves on 
surface. However, different from traditional periodic and sohtary waves with smooth 
crest, the horizontal bottom velocity of the peaked solitary waves is always larger than 
that on free surface. Certainly, due to the viscosity of fluid, the horizontal velocity of 
all water waves must be zero at bottom, so that such kind of the peaked solitary waves 
might not exist exactly in its theoretical form reported in this article, since there exists 
a thin viscous boundary layer near the bottom. Thus, such kind of solitary waves with 
a peaked crest and larger bottom velocity might be more difficult to create not only 
in laboratory experiments but also in nature. This might be the reason why such kind 
of peaked solitary surface waves have never been reported and observed. However, if 
such kind of peaked solitary waves with larger horizontal velocity near bottom could 
indeed exist in nature, even though they might be not exactly in the same form, it 
would be quite dangerous to submarines, platforms and equipments in underwater 
engineering. 

Note that the peaked solitary waves found in this article are obtained under the 
assumptions that the fluid is inviscid and incompressible, the flow is irrotational in 
the domain x > (it implies that the flow is not necessarily irrotational at x = 0), 
the surface tension is neglected and the wave elevation has a symmetry. Although 
the same assumptions are widely used for the traditional periodic and solitary waves 
with smooth crest, their physical reasonableness for the peaked solitary waves should 
be reconsidered carefully in future. Especially, the surface tension might have an 
important influence on the crest of the peaked solitary wave. Besides, the viscosity of 
the fluid also has an important influence on the horizontal velocity at bottom and the 
discontinuous vertical velocity under crest (i.e. at a; = 0). These influence might be 
local and might not essentially change the properties of the peaked solitary waves. All 
of these are interesting and should be investigated in details. Without doubts, further 
theoretical, numerical and experimental studies, and especially practical/experimental 
observations about this new type of solitary surface waves with peaked crest and 
many unusual characteristics, are needed in future. All of these might deepen our 
understandings and enrich our knowledge about solitary waves. 

Finally, it should be emphasized that, both of the smooth base-functions (I16p and 
the so-called evanescent-mode base-functions fl2B]) are familiar and can be found in 
the textbook of water waves [1A\. Both of them satisfy the linearized wave equations. 
Especially, waves expressed by each kind of these base-functions admit the disconti- 
nuity of the velocity. Besides, the discontinuity of wave elevation exists widely not 
only in theory, such as Stokes limiting gravity wave with a corner at crest and the 
peaked (see 0) and cusped [12] solitary waves of the CH equation ([T]), but also in 
natural phenomena, such as dam break |52| - [55] in hydrodynamics, shock waves in 
aerodynamics, and so on. Therefore, the peaked solitary waves reported in this ar- 
ticle are consistent with the traditional waves with smooth crest. Besides, it also 
provides a reasonableness of the existence of the published peaked or cusped solitary 
waves IS1II2]- However, the peaked solitary waves reported in this article have many 
unique, unusual characteristics different from the traditional waves with smooth crest, 
which clearly indicate their novelty. 
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Indeed, the discontinuity and/or singularity are difficult to handle by traditional 
methods. But, they should not be evaded easily, since they might open some new 
fields of research, and greatly enrich and deepen our understandings about the real 
world. 
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Appendix 

Detailed derivation of formulas (I74[) - (1791) 



In this appendix, we explicitly give the formulas of all terms in ( FM|) - (1791) . 
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Then, on the unknown free surface z = t]{x; q), we have using (!97|) that 
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+00 

dx^ V dz 



where 



7. 



ji.O 
n,m 



V 92;^ 



m=0 

+00 



m=0 



rn,0 



E(^5 + l)(^ + 2)C''^Vs,,n, m>l. 



s=l 



Then, on the unknown free surface z = rj{x; g), it holds using fllU3l) that 



+00 



+00 



<l>(x,C;g) = E'^"(^'0 9" = E^" 



n=0 

+00 +00 



n=0 



71=0 m=0 

+00 



+00 

E /^o '"(^) <r 

.m=0 

+00 

s=0 



m=0 



s—m,m I 



n=0 



where 



.0(3;) = E^' 



n—m,m 




45 



Similarly, we have 
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Then, on the unknown free surface z = r]{x;q), it holds using flll4p and (I115p 
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and similarly 
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Substituting f ll34p . ( I128P into ( 158|) and equating the like-power of g, we have the 
boundary condition: 
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Substituting the above expression into fll36p gives the boundary condition on z = 0: 
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Substituting the series fl70|) . flll4p and fll20p into fl59l) . equating the like-power of 
g, we have on z = that 
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